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Cluster structures of excited 10 Be states are investigated with a hybrid model of dineutron conden- 
sate wave functions and 6 He+a cluster wave functions. Two kinds of cluster states are theoretically 
suggested a few MeV above the a+a+n+n threshold energy. They have quite distinct cluster struc- 
ture; ones have gas-like structures of a+a+dineutron, and the others have 6 He+a with extremely 
extended an a cluster. Although these cluster states have not been confirmed experimentally yet, 
theoretically suggested properties of these states such as monopole transition strengths and a-decay 
widths may be helpful for possible experimental observation. 



I. INTRODUCTION 

Nuclei have a variety of unique features. A cluster 
structure is one of them, which means that one or a few 
subsystems composed of some nucleons are formed in nu- 
clei. Nowadays cluster structures have been investigated 
eagerly both theoretically and experimentally, and many 
efforts have been made to understand their origin and 
features. The most representative cluster is an a particle 
which is very stable and formed in a number of nuclear 
systems. 

Be isotopes are the typical systems where a clusters de- 
velop. The structures of Be nuclei heavier than 8 Be can 
be described so well by a 2a core and valence neutrons 
around it fJ-TTol] - 10 Be is the lightest even-even neutron- 
rich Be isotope and shows diverse structures composed of 
2a plus two valence neutrons, such as molecular orbital 
structures. Such structures have been investigated by 
many theoretical methods such as molecular orbital mod- 
els p], H, d, [U , cluster models [H, Q and antisymmetrized 
molecular dynamics [3,[I3] and so on. As for its + states, 
the and 0j states have been confirmed experimentally 
and their structures are well-described theoretically. Re- 
cently, in addition to those two + states, a possibility of 
the existence of a new + state was suggested experimen- 
tally on the analogy of 10 B According to the report, 
the + state should be located about 3 MeV higher above 
the 6 He+a threshold and might have a large a spectro- 
scopic factor. 

In relation to cluster structures in nuclear physics, we 
are interested in dineutron correlation in finite nuclei. A 
dineutron is composed of two neutrons coupled to spin- 
singlet which are not bound in free space. It is theo- 
retically suggested that dineutron correlation is rather 
enhanced in the low-density region in the nuclear matter 
fl2l [HI ] and at the surface of finite nuclei, especially the 
neutron-halo and -skin region of neutron-rich nuclei 
[l7j . So, even if a dineutron is a fragile subsystem unlike 
an a, it is expected that dineutron clusters are enhanced 
in nuclei, especially at the surface of neutron-rich nuclei. 
It is also a challenging problem to search for developed 
dineutron structures in excited states as one suggested 
in 8 He [l8|. In order to investigate dineutron correla- 
tion in finite nuclei, we had constructed a model, which 



we name the dineutron condensate (DC) wave function, 
and applied to 10 Be In that work, we suggested a 
new + state composed of a 2a core plus a dineutron 
above the 0^" state, though we could not come to a con- 
clusion about the detail of that state. The possibility 
of the 2a core plus a dineutron structure has also been 
discussed in Ref. [20|. Then, a question is whether the 
theoretically suggested 2a core plus a dineutron state in 
Ref. [l9[ can be assigned to the experimentally suggested 
state in Ref. [ll[ or they are individual states above the 
threshold energy. 

In the present work, in order to investigate cluster 
structures of 10 Be in detail from the viewpoint of struc- 
tures of 6 He+a and a+a+dineutron, we superpose two 
kinds of cluster wave functions, 6 He+a cluster wave func- 
tions and a+a+2n DC wave functions. The developed 
6 He+a cluster states are described in detail by the 6 Hc-cv 
cluster wave functions and structures containing a dineu- 
tron cluster are expected to be described by the DC wave 
functions. Through such a close investigation, we will 
conclude that the cluster states mentioned above would 
exist as distinct ones. In this paper, we mainly discuss on 
the characteristic structures of those cluster states which 
have not been established experimentally yet. We show 
that the ones have almost independent a+a+dineutron 
clusters as suggested in our previous work, and that the 
others are the developed 6 He+a cluster states, the + 
state of which, we believe, corresponds to the state sug- 
gested in the experimental report of Ref. [TTI | . We also 
discuss such properties of these states as monopole transi- 
tion strengths and a-decay widths, which would be help- 
ful information for experiments to search for these states. 

The outline of this paper is as follows. In Sec. lU we 
present our framework in this work briefly. In Scc. lIIIl we 
show the obtained results and discuss features of the clus- 
ter states. In Sec. IIVI we summarize the present work. 



II. FRAMEWORK 

In the present work, we superpose two kinds of cluster 
wave functions in order to describe 10 Be states, 6 Hc+a 
cluster wave functions used in Ref. 2l| and 10 Be DC wave 
functions used in Ref. [l^]. The details of the formalism 
for each model are explained in these references and we 
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do not enter the details here. 



A. The He+Q cluster wave function 

In 6 He+a cluster wave functions, |$ c (d c )), the 6 He 
and a clusters are separated by a distance d c . The 
spatial parts of the single-particle wave functions in a 
and 6 He are described by the harmonic oscillator (H.O.) 
shell model wave functions whose centers of mass are lo- 
cated at (0,0,3ci c /5) and (0, 0, — 2d c /5) respectively and 
widths are the same value of v = l/(26 2 ) = 0.235 fm -2 . 
In an a cluster, four nuclcons are in the s-shell closed 
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The 6 He cluster is composed of an a 



core and two valence neutrons in Op-orbits around the a 
core. We consider all the Op-shell configurations for two 
valence neutrons and mix those configurations to describe 
the 6 He cluster. We describe the Op-configurations with 
the shifted Gaussians from the a core as in Ref. [2l[ . An 
isolate 6 He described in such the simple configurations 
of the H.O. shell model does not necessarily agree with 
the realistic one qualitatively. For example, its binding 
energy is higher than the threshold energy to a+n+n 
and it has no neutron-halo structure. In spite of such a 
simplified 6 Hc, the structure of the well-known states of 
10 Be can be described well with the 6 He+a cluster wave 
functions as discussed in Ref. 1211 . 



B. The 10 Be DC wave function 

In addition to the 6 He+a cluster wave functions, we 
superpose the DC wave functions containing dineutron 
correlation explicitly. The 10 Bc DC wave function is com- 
posed of the determinant of the eight single-particle wave 
functions in the 2a core and the two ones in a dineutron 
which contains two neutrons coupled to spin-singlet. Its 
spatial form is 
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In Eq. ([!]), the single-particle wave functions in an a clus- 
ter, ip a i (i = 1,2), are the same as used in 6 Hc+a cluster 
wave functions, whose centers are at (0,0, ±d n /2) and 
widths are 0.235 fm~ 2 . The wave functions of two neu- 
trons in the dineutron have the common Gaussian width, 
b n , and center, Y n . The center Y n is integrated out 
multiplied by the Gaussian weight. The width of the 
Gaussian weight is B n and the dineutron is distributed 
around the 2a core in the spherical 5-orbit on the whole 



space. In order to make clear the relative motion be- 
tween two neutrons and between the dineutron and the 
core, we perform the l^„-integral and rewrite two neu- 
tron wave functions to the relative and center of mass 
wave functions of two neutrons as follows. 
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Here we define a new parameter 1 = B 2 + 6 2 . The 
parameter b n is the width of the relative wave function 
(Eq. (HI), and (3 stands for the one of the dineutron cen- 
ter of mass wave function (Eq. ([S])) so that we call the 
parameters b n and /3 as the dineutron size and the dineu- 
tron extension from the core respectively. In DC wave 
functions, we take various values for the parameters b n 
and (3 and superpose the DC wave functions to describe 
the changing effect of the dineutron structure as well as 
the one of the dineutron size. It should be noted that 
the Gaussian width for two neutrons in a dineutron, b n , 
is generally different from the one of the single-particle 
wave functions in the core, 6, and as a consequence, it 
is difficult to separate the center of mass motion exactly. 
So, in calculating energy expectation values, we approxi- 
mately treat the center of mass motion effect by extract- 
ing the expectation value of center of mass kinetic energy 
from the total energy. On the other hand, in calculating 
root mean square radii and so on, we neglect the effects 
of the center of mass motion for simplicity. 



C. The 10 Be wave function in the present work 



By superposing the 6 He+o; cluster wave functions and 
10 Be DC wave functions and projecting each wave func- 
tion to the parity and angular momentum eigenstates, 



we describe °Be(J 7r ) states, j'S 
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Here, the indices i, j are the abbreviations of a set of the 
parameters, that is, i = {d c } and j = {d n ,b n , f3}, and 
the coefficients cxi,Kj arc determined by diagonalizing 
the Hamiltonian. In the present work, we superpose the 
6 He+a cluster wave functions of d c = 1,2, •••8 fm to 
describe the structure of 6 He plus a rather developed a 
cluster. We also superpose the DC wave functions of five 
b n values, j3 = 2, 3, ■ • • 9 fm and d n = 1, 2, ■ ■ ■ 6 fm. The 
set of b n values are chosen for each /3 in the same way 
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as in Rcf. [19|. Wc superpose DC wave functions with 
rather large (3 (up to 9 fm) in order to describe the state 
of interest which contains a weakly interacting dineutron 
with respect to two as. Such constraints correspond to 
a bound-state approximation with respect to a+a+n+n 
and 6 He+a. 

In the present work, we investigate only the even-parity 
states because the DC wave functions of 2a plus a dineu- 
tron contain only even-parity components. 



III. RESULTS 



A. The effective Hamiltonian 



In this work, we use the effective Hamiltonian as 



H = T-T G + V ccnt + V LS + V c 



'Coul, 



(7) 



where T and Tq are the total and center of mass kinetic 
energy. We used the Volkov No. 2 force [22[ as the central 
force, V^nt, and the spin-orbit part in the G3RS force [23| 
as the spin-orbit force, Vls- The Coulomb force, Vcoub 
is approximated by seven-range Gaussians. The adop ted 
parameters are the same as those used in Refs. |lOl liH 
l2l| . that is, in the central force m = 0.60, b = h = 0.125 
and in the spin-orbit force the strength v\ = —v-2 = 1600 
MeV. We also use modified parameters below to see the 
dependence of the level structure of the states of interest 
on the parameter choice. 



B. Energy spectra 

The energy spectra of the even-parity states obtained 
by the present calculation compared with experimental 
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FIG. 1: Energy spectra of experiments and the present cal- 
culation on the basis of the ground state energy. The exper- 
imental and calculated threshold energy to a+a+n+n and 
to 6 He+a are shown with dotted and dashed-dotted lines re- 
spectively. The assigned spin-parity, J w , are labeled to each 
state and the orders only in our result. 



one are shown in FIG. [I] Experimentally established + 
states are the ground + state and the first excited + 
state at E x = 6.2 MeV, which exist 8.4 MeV and 2.2 
MeV below the a+a+n+n threshold, respectively The 
0^2 states obtained in the present results correspond to 
these two + states. In the calculation, the energy posi- 
tions of the Of and 0^ states relative to the a+a+n+n 
threshold energy are overestimated compared to the ex- 
perimental values. This discrepancy mainly originates in 
the fact that in the present calculation the energy of 6 He 
is overestimated and the 6 He+a threshold is much higher 
than the a+a+n+n threshold in contradiction to exper- 
imental data. Such a disagreement is also seen in other 
calculations with similar effective interactions 0, [E [HI ■ 
In spite of the failure in the quantitative reproduction 
of the energy position, the energy levels and structure 
properties of the states in 10 Bc have been reasonably re- 
produced by the pioneering works. As shown later, the 
disagreement can be improved by slight modification of 
interaction parameters, which gives qualitatively similar 
results in the present calculation. 

In addition to the 0^~ 2 states, two + states above the 
Of state are obtained in the present calculation, which 
we tentatively label 0^ 4 states here. Each of them forms 
the rotational band with the members of + , 2 + and 4 + 
states. In the 0j band, the angular momentum J is pro- 
duced by the rotation of 2a relative motion, while in the 
O4" one, it is mainly contributed by the angular momen- 
tum of 6 He and a relative motion. The members have 
the same structure with those of their band head states, 
so we focus on the O3 4 states in this paper. As shown 
closely in the next subsection, the J state is mainly com- 
posed of developed two as and one dineutron cluster and 
the O4" state has mainly 6 He(0 + ) with an extremely de- 
veloped a. Since the 0^" and 0^ states are expected to be 
quasi-bound states, they should be coupled to continuum 
states of a+a+n+n and 6 He+a respectively. In fact, al- 
though these states are calculated within a bound-state 
approximation as mentioned in the previous section, they 
arc still coupled to the neighboring continuum states to 
some extent and the continuum states have some frag- 
ments of the 0^4 states. In the present paper, we consider 
the states having the most remarkable cluster structures 
which wc show in the next subsection as each represent of 
0^4 without a width for simplicity. In addition to consid- 
ering the significant cluster components in each state, we 
also check the identification of those quasi-bound state 
by using the pscudo potential and separating from con- 
tinuum states as explained in Appendix. 

Here it should be noticed that the 6 He+a threshold is 
much higher than the a+a+n+n threshold in the present 
calculation, as is different from the experimental fact that 
they almost degenerate. The failure originates in that 
it is generally difficult to reproduce systematically the 
binding energies of subsystems by effective two-body nu- 
clear interactions. In this work, we focus on the 6 Hc+a 
cluster states and the 2a and a dineutron cluster states, 
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the qualitative features of the structure of those states 
are almost unchanged. Below, we discuss the results ob- 
tained by the original interaction set. 



We investigate the detailed structures of the 0f 4 states, 
which we suggest are characterized by developed two as 
and one dineutron clusters, and extremely separate 6 He 
and a clusters, respectively. 



default m=0.57 + b=h=0.0 + v LS =2500 



The 0j state of two as plus one dineutron 



FIG. 2: Energy spectra calculated with different parameter 
sets. Each set is explained in the text. The correspondent 
states are linked with dotted lines. And the threshold energy 
to a+a+n+n, a constant of —55.23 MeV, and the ones to 
6 He+a, which depends on the parameters, are shown with a 
thin solid and dash-dotted lines respectively. 



which might be sensitive to the corresponding thresh- 
old energies, and therefore, the relative position of these 
threshold energies can be important to conclude the level 
structure of the 0^" 4 states. To ensure the interaction de- 
pendence of the results, we calculate the + energies by 
using four kinds of different sets of parameters. The en- 
ergy spectra with the different parameter sets are plotted 
in FIG. [21 In the first column which is labeled as "de- 
fault", we plot the spectrum calculated with the orig- 
inal parameter set of m = 0.60 and b = h = 0.125 
in the Volkov No. 2 force and Vi = -v 2 = 1600 MeV 
as the strength in the spin-orbit force, as mentioned in 
Sec. MI Al In the second labeled "m=0.57", we calcu- 
late by changing only the strength of Majorana exchange 
term to m = 0.57. The spectrum with b = h = 0.0 for the 
strengths of Bartlett and Heisenberg terms and m = 0.57 
is shown in the third column (+ b=h=0.0). Besides the 
set (6 = h = 0, m = 0.57) in the third, we modify the 
strength of the spin-orbit force to 2500 MeV and show 
the calculated spectrum in the last column (+v L s=2500). 
By such a series of the change in the parameters, the 
binding energy of 6 He decreases but the a energy is un- 
changed, so that the 6 He+a threshold position relative 
to the a+a+n+n threshold is improved substantially. It 
is found that the level ordering is qualitatively similar for 
these four parameter choices. In particular, the energy 
positions of the 2 and Of states relative to the 6 Hc+a 
energy is not sensitive to the parameter choice. On the 
other hand, the Of state, which have a quite distinct 
component from those of the 0f 4 states, shows a little 
different dependence. So it would be possible that, as 
shown in the fourth column, the 0^" 4 states are located 
at closer energy positions, that is, a few MeV above the 
threshold to a+a+n+n or to 6 He+a. 

In spite of the interaction dependence of the energies, 



At first, we discuss the dineutron component in the 0g 
state. We show the overlaps with DC wave functions, 
A/jn; a s the components of a dineutron cluster, 

M 2n (b n ,(3,d n ) = |<< c (&„,/Mn)|* 0+ >| 2 , (8) 

plotted on the j3-b n plane for fixed values of d n in FIG. [3] 
A component with a small b n corresponds to the one of 
a compact dineutron, and a component with large /3 to 
the one of two neutrons distributed far from the 2a core, 
and vice versa. The significant broad peak is seen at the 
region of 3 < (3 < 6 and 1.5 < b„ < 3. These major peaks 
correspond to the component of a rather compact dineu- 
tron far from the core. The significant amplitudes of the 
peaks existing in the d n = 3, 4, 5 fm cases indicate that 
the a-a distance fluctuates in the 0^ state. The minor 
peak around (P,b n ) ~ (1.5,1.5) corresponds to the one 
of two neutrons distributed in the middle of two distant 
as. Since such the ((3,b n ) ~ (1.5, 1.5) component is the 
major component in the Of state as shown in Ref. (l9j |. 
the amplitudes around the minor peak in the 3" state are 
thought to come from the orthogonality to the Of . These 
behaviors suggest that two a and one dineutron clusters 
are weakly interacting to each other and spread over a 
broad space, that is, behave as a gas-like state consisting 
of two a and one dineutron clusters. 

We should offer a comment on a difference between 
the behavior of a dineutron in the 0^ and Of states. As 
discussed in Ref. the Of state has components of 
a dineutron whose peak is at (/3,b n ) ~ (2.0, 1.5), which 
corresponds to a very compact dineutron just at the nu- 
clear surface. On the other hand, the 0^ state has a main 
peak at (f3, b n ) ~ (4.0, 2.3) in FIG. [3fb), where a dineu- 
tron has a larger size and a broader expansion from the 
2a core than those in the Of state. We suppose it can 
be said that a dineutron is so fragile that its size is very 
changeable according to the very nuclear structure. 

As mentioned above, the 10 Bc(0 3 f ) shows the behaviors 
of weakly interacting two a and one dineutron clusters, 
which can be associated with the 3a-cluster gas state 
suggested in the Hoyle state, 12 C(0^~), by replacing one 
of three as with a dineutron cluster. However, it is ex- 
pected that the dineutron cluster is not as rigid as an 
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a because, as well known, two neutrons are not bound 
in free space and even in nuclear media a dineutron is 
naturally expected to be a rather weakly (semi)bound 
subsystem. In particular, the dineutron size may enlarge 
as the dineutron goes away from the core. Indeed, the 
representative b n value at the peaks, ~ 2.3 fm, is larger 
than the size of an a, ~ 1.5 fm. Furthermore, it is seen 
in FIG. [3] that the tail part of the main peak spreads to 
large (3 region. It suggests that two neutrons would ex- 
pand broader and the lo Be(0^~) may have a broad width 
inevitably. 



2. The Oj state of 6 He plus a 

Next, we discuss 6 Hc+a cluster structure of the O4" 
state. To see how far an a cluster is separated spatially 
from 6 He, we analyze the overlap with 6 He+a cluster 
wave functions, 

N a (d a ) = |<P 0+ $a H e(o+) ® ® a (d a )\* 0+ )\ 2 . (9) 

Here the wave function $6Hc,a is fixed at the cluster dis- 
tance d a and the total wave function is normalized and 
projected to the + state. The 6 He ground state wave 
function ^ > 6 He (Q+) is approximated to be the lowest state 
calculated in the H.O. p-shell configurations by using the 



present effective interactions. When d a is large, it corre- 
sponds to the components of a developed a cluster plus 
the ground state of 6 He(0 + ). 

The overlaps of the 0^ 4 states are shown in FIG. [4j It 
is found that the state has striking amplitude in the 
large d a region, whose peak is at about 7 fm. It indi- 




d a (fm) 



FIG. 4: The overlaps of the O3 4 states with He+a cluster 
wave functions, M a , which are plotted with a dotted and solid 
line respectively. The horizontal axis is the distance between 
6 He and a. 
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cates an extremely developed a cluster in the Of state, 
and that is consistent with the experimental indication 
in Rcf. [Til] - Such a large distance as 7 fm is almost 
out of the inter-nucleus potential range. In back of the 
mechanism of the formation of the Of state, a rather low 
Coulomb barrier exists around 8 fm from 6 He. The a 
cluster is confined in such a wide region and constitutes 
the so-called "a-halo" referred in Ref. 24 1. This state 
should be measured by a-decay due to such a developed 
a cluster. 

In contrast to the case of the Of state, the Of state 
shows not significant but only tiny amplitude in the d a < 
6 fm region. This state has developed two as and one 
dineutron. In such the state with the weak correlation 
between as and a dineutron, a 6 He cluster cannot be 
formed. It can be a reason why the component of a single 
separated a from the 6 He cluster is suppressed in the Of 
state. 

More detailed analysis of inter-cluster motion between 
6 He and a clusters will be given in the next subsection. 



D. Discussion on experimental observables 



section. In the 6 He(0f ) state, the (OP3/2) 2 configuration 
of the two valence neutrons is dominant (~ 99%) in the 
present calculation. 

The calculated RWA of each + state is shown in 
FIG. In FIG. EJa), it is seen that the Of state has 
very large amplitude in the large r region as expected. 
Compared to the RWAs of the Of and Of states, the 
node number is the minimum in the Of and it increase 
one by one in the Of state and in the Of state. It sug- 
gests that the Of state might be interpreted as the higher 
nodal state of the Of state caused by the excitation of 
the inter-cluster motion. Incidentally, its amplitude is 
rapidly decreased over 8 fm, for the present model space 
is truncated by d c = 1 — 8 fm 6 He+a cluster wave func- 
tions. As mentioned in Sec. Ill C[ the present 10 Be wave 
functions are obtained within the bound state approxi- 
mation against 6 He+a continuum states. The Of state 
actually should couple to the continuum states further 
and it might have larger amplitude in the large r region. 
To see the effect of further mixture of 6 He+a continuum 
states, we compare the RWAs of the Of state calculated 
by superposing cluster wave functions with d c = 1 — 8 
fm and d c = 1 — 10 fm in FIG. [Hfb). In the case of 



As mentioned above, these Of 4 states have distinct 
structures from each other, and it is expected that they 
can be distinguished in experiments by measuring char- 
acteristic quantities. Here we discuss the a-decay widths 
and the monopole transition strengths for the Of 4 states, 
and show that the a-decay would be useful to confirm the 
Of state and that the monopole transition is promising to 
observe the Of state in inelastic scattering experiments. 



1. The a-decay width 

Since an a cluster is developed in the exotic cluster 
states of interest, especially in the Of state, they should 
be confirmed by surveying a decay to 6 He and a. There- 
fore we evaluate the partial a-decay widths by using the 
reduced width amplitudes (RWAs). The RWAs are calcu- 
lated following an approximate method in Rcf. Q where 
the eigenvalues of the norm kernel are calculated numer- 
ically. 

The RWA between 6 He and a, 3^, and the a spectro- 
scopic factor, S a , are defined as follows. 



s aL = / dr r z y L ( r y, 



(10) 

(11) 



where a is the channel radius and Ylq(t) is the spher- 
ical harmonic function. </>o( 6 He, a) is the internal wave 
function of 6 He(0f ) and a cluster and $z, is a 10 Be wave 
function. Here the 6 He(0f ) wave function is the inter- 
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FIG. 5: (a)The RWAs of 0+ 
axis is the distance between 



nal wave function of $ 6Hc (o+) explained in the previous wave functions are up to 8 fm or 10 fm. 
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TABLE I: The a-decay widths, F a , and the a spectroscopic factors, S a , of the 0^" 4 states in the cases where the resonance 
energy, E r , is 2,3,4 MeV. In the calculations of the Q-decay widths, the channel radius, a, is chosen to 5, 8 fm for 0^4 respectively. 
a spectroscopic factors do not depend on the resonance energy so that we write the values only on E r = 2.0 MeV lines. 







d c = 1 - 8 fm 




d c = l- 


- 10 fm 




E r (MeV) 


T tt (keV) 


Sot 


T Q (keV) 


Sot 




2.0 


2.7 x 10 1 


0.03 


2.8 x 10 1 


0.03 




3.0 


3.8 x 10 1 




4.0 x 10 1 






4.0 


4.7 x 10 1 




4.9 x 10 1 




oi 


2.0 


3.5 x 10 2 


0.21 


4.2 x 10 2 


0.28 




3.0 


4.7 x 10 2 




5.6 x 10 2 






4.0 


5.6 x 10 2 




6.6 x 10 2 





d c = 1 — 10 fm, the amplitude in the large r region is 
actually enhanced because of coupling with the contin- 
uum states. However, the amplitudes around the peak at 
r = 7 ~ 8 fm does not change so much, indicating that 
the present bound state is reasonable. 

Using the RWA derived as above, the partial a-decay 
width, r„r, is estimated by a low-energy limit approxi- 
mation 125! as follows. 



P L {ka 



T aL = 2P L (ka)j 2 aL (a), 
ka 



Fl(ka) + Gl(kaY 



(12) 
(13) 

(14) 



where [i is the reduced mass of 6 He and a, and k 
is the wave number of the resonance energy E r (k = 
y / 2^iE r /h 2 ). The partial a decay width, T a L, is the prod- 
uct of the reduced width, 7 q l, and the penetrability, Pl- 
In Pl, Fl and Gl are the regular and irregular Coulomb 
wave functions respectively. 

The a-decay widths and the a spectroscopic factors of 
the 0^ 4 states are shown in TABLE [TJ In the calculations 



of the a-decay widths, we use the channel radius, a = 
5, 8 fm for the 0^" and O4" states respectively, which we 
choose to slightly larger values than those of the main 
peak positions of the RWAs shown in FIG. [5^a). For 
the resonance energy, we use E r = 2,3,4 MeV for both 
states, which are near the energy suggested in Ref. [ll[ as 
the resonance energy above the threshold to 6 Hc+a for 
the possible 6 He+a resonance state. The 0^ state has 
the larger decay width and a spectroscopic factor than 
those of O3 as expected. When the 6 He+a cluster wave 
functions up to d c = 10 fm are superposed, the width 
and a spectroscopic factor are enhanced largely in the O4" 
state because this state is very sensitive to the component 
of separated 6 He and a, while calculated values hardly 
change in the O3" state, which contains a little component 
of 6 He+a. It indicates that the a-decay width of the 0^ 
state shown in TABLE. UJ would still become larger by 
considering the coupling to continuum states fully. It is 
noted that the approximation, Eq. (|12p , is appropriate 



for narrow resonances but may be inappropriate for such 
broad ones as the Ot states. To estimate the broad a- 
decay width of the O4 the state more correctly, we use the 
formula which is the origin of the approximated Eq. (|12|) 
referred in Ref. 1261. 



Gl 
ka 



1 



Sl-Pl 



Fl 



(15) 



where S and 4> are the usual shift function and the hard- 
sphere scattering phase shift, and the dots signify the 
energy derivative. Applying this formula to the O4 state 
formed by superposing d c = 1 — 10 fm 6 He+a cluster wave 
functions, we obtain T a — 8.2 x 10 2 keV for E r = 3.0 
MeV and a = 8.0 fm, which is about 50 % larger than 
the r a evaluated with Eq. (ri"2|) . 

Besides, although we show here the results of the rep- 
resentative states having the major component of each 
state, part of the components of the 0^~ 4 states are frag- 
mented in neighboring continuum states. Therefore, the 
actual a-decay width which is the sum of those fragments 
would be somewhat larger than our calculated values. It 
should be also noticed that we do not consider a one- and 
two-neutron decays, and therefore, the total decay width 
should be still larger. In order to calculate the more ac- 
curate values for the widths, it is essential to perform a 
more sophisticated treatment of a resonance with a broad 
width considering all the decay channels, but it is beyond 
the present study. We just emphasize that the O4 state 
certainly has the large a-decay width and it can be con- 
firmed by surveying a-decay. 



2. The monopole transition strength 

The monopole transition is expected to be a measure 
of the extension of clusters [13, HH- Two factors are 
essential for a large monopole transition strength. The 
first is that states before and/or after transition have a 
extended structure because the monopole transition op- 
erator is composed of r 2 (Eq. (flB|) ) so that it reflects the 
degree of a radial excitation. The second is that the con- 
figuration docs not so differ before and after transition, 
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TABLE II: The total, proton and neutron monopole transition strength from the state to the ground state, Ml' p ' n . In the 
parentheses, the ratios to the first-order EWSR, 5 , 1 ' p ' n , are shown by %. 





Mt 




MS 




MS 






5.9 


( 5.3 % ) 


2.1 


( 1.8 % ) 


3.8 


( 3.5 % ) 




8.6 


( 13.6 % ) 


1.5 


( 1.1 % ) 


7.1 


( 14.7 % ) 


ot 


2.9 


( 2.2 % ) 


2.4 


( 4.1 % ) 


0.5 


( 8.8 xl(T 2 % ) 



for example, it can be large in the case of the transition 
from an excited cluster state to a ground state which in- 
volves some cluster components. So it is expected that a 
remarkable monopole transition strength can be observed 
in 10 Be where the ground and excited states have striking 
cluster structures. 

The total monopole transition strength between the 
excited 0+ states and the ground state, Ml, is defined as 

Mt = (0t\J2rf\0t), (16) 

i 

where the index i runs over all of the nucleons. The pro- 
ton and neutron ones, M^ n , are denned in the same way, 
though the index i runs over only protons or neutrons. 
Here these calculations are performed by superposing DC 
wave functions whose /3 is reduced to 2 — 6 fm and 6 Hc+a 
cluster wave functions with d c = 1 — 8 fm. Both of the 
0^4 states are spatially extended quasi-bound states so 
that they tend to be mixed with continuum states of 
a+a+n+n and 6 He+a respectively. The more contin- 
uum states are coupled, the more strength increases due 
to their contributions. Therefore, with the inclusion of 
a larger model space and more coupling with continuum 
state, calculated monopole transition strengths become 
very uncertainty So in the calculations of monopole 
transition strengths, we adopt the further restricted DC 
wave functions. We have confirmed that their qualitative 
structures do not differ from those of the original ones. 
As the criterion of a degree of the strength, we also cal- 
culate the ratios to the first-order energy-weighted sum 
rule (EWSR), 

Sl = £(^-^)l<0j|£r?|0+)| a , (17) 

v i 

= l:Kl£«>> as) 

i 

where E v — Eq is the excited energy of the ^th + state. 
The second line, Eq. (fl"8|) . holds in the case that the 
Hamiltonian contains only momentum-independent in- 
teractions. The definitions of the superscripts t,p,n are 
the same as those in the monopole transition strength. 
The calculated values of the strengths between the 
ground state and the 0^ 3 4 states and their ratios to the 
EWSR are shown in TABLE M The j state has a re- 
markable strength in the neutron one and it exhausts 14.7 
% of the EWSR. It reflects the remarkable gas-like struc- 
ture of the 0^" state shown in FIG.[3J This state has two a 



and one dineutron clusters which are weakly interacting 
with each other so that the neutrons are spatially devel- 
oped largely. Moreover, the ground state contains two 
a clusters and remarkable dineutron correlation at the 
surface as shown in Ref. As a result, its monopole 
transition strength is highly enhanced. In contrast to the 
$ , the state has just a little strength of a few % of the 
EWSR. It is because the O4" has the developed a cluster 
structure which is so extended that the overlap with the 
ground state is too small. Indeed, the RWA of the O4 has 
the second higher nodal structure of that of the state 
and shows different behaviors (FIG. Ufa)). Thus the dis- 
tinct structures between the 0^ and 0^ states cause the 
clear difference in their monopole transition strengths, 
which would be a good probe to identify the and O4 
states. It is indeed meaningful that not the absolute val- 
ues but the relative magnitude since we superpose only 
restricted wave functions in the present calculation of the 
monopole transition strengths. The present result sug- 
gests that the 0^" state has the more striking monopole 
transition strength than those of the other + states, so 
that we expect that the 0$ state can be confirmed ex- 
perimentally in the inelastic scattering from the ground 
state. 



IV. SUMMARY 

We have investigated the exotic cluster structures of 
excited states of 10 Bc by using the 6 He+o; cluster wave 
functions plus the DC wave functions. In the present 
study, we suggest theoretically two kinds of novel cluster 
states above the a+a+n+n threshold energy, which have 
not been confirmed experimentally yet. These states have 
remarkable cluster structures and construct the K* = + 
rotational bands which show quite distinct characteris- 
tics. Ones of them contain developed two a and one 
dineutron clusters. They are weakly interacting to each 
other, and the + state in the band is regarded as a 
gas-like structure in association with the Hoylc state, 
12 C(0j). In addition to such states, we also found states 
composed of 6 He and a developed a cluster. Because of 
a large a spectroscopic factor, the + state is likely to 
correspond to the one suggested experimentally on the 
analogy of 10 B. In this state, the a cluster is extremely 
extended from 6 He so that one could even call it a-halo. 

Since these new states have distinct cluster structures, 
they should be distinguished with the different observ- 
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ables. We have suggested that the + state including 
extended two as and one dineutron can have a large 
monopole transition strength from the ground state due 
to the gas-like extended a+a+dineutron structure. On 
the other hand, the 6 He+a state has an extremely devel- 
oped a cluster so that it can be confirmed by measuring 
the decay to 6 He and a. We have made sure that the 
6 He+a state has a much larger a-decay width than the 
a+a+dineutron state. 

In the present work, we have mainly focused on the 
structures of two novel cluster + states of 10 Be. These 
states are obtained within the bound state approxima- 
tion and are suggested to be resonance states above the 
a+a+n+n and 6 He+a threshold energies. However, in 
order to predict precise resonance energies and widths of 
these states, further improvements of the calculation are 
required. For instance, we need to treat resonance behav- 
iors beyond the bound state approximations by taking 
into account all the decay channel which are neglected 
here and coupling with continuum states. Moreover, 
the threshold energies, in particular, those of a+a+n+n 
and 6 He+a channels should be reproduced by improv- 
ing wave functions and adjusting effective interactions 
carefully. The purpose of this work is to suggest exotic 
cluster states and the points mentioned here are tasks for 
the future. 

Dineutron correlation is usually discussed in study of 
ground state properties. What we have shown in this pa- 
per is that it can be important not only in ground states 
but also in excited states. As a future work, we will ap- 
ply our method systematically to the nuclei neighboring 
10 Be, that is, 9 Li and 8 He, and search for states contain- 
ing one or a few dineutron clusters and investigate the 
tendency of dineutron cluster formation in those nuclei. 



and add it to the original Hamiltonian, H, (Eq. ([7|), 



H(6) = H + S xV 
V = vq exp 



i<j 



a 



(A.l) 
(A.2) 



where as the parameters v = —100 MeV and a = 1.0 fm 
are chosen. S is the parameter to control the strength of 
the pseudo potential. When 8 — 0, the modified Hamil- 
tonian, H(S), defined in Eq. (jA.ljl is equal to the original 
one, H. When S is increased, since the short-range at- 
traction between nucleons becomes artificially larger, res- 
onance states which have larger components in a finite re- 
gion should gain relatively more energy than continuum 
states. With an enough large 6, the resonance states 
come down below continuum states to become bound 
states. As a result, the states of interest can be sepa- 
rated out of continuum states energetically. 

By diagonalizing the Hamiltonian matrix calculated 
with H(S), energy levels are obtained. The energies of 
+ states calculated as below are plotted in FIG. [6l 

AE(6) = ^ 0+ (S)\H(6)\^ + (S))-^t(S)\H(S)\^f a (S)). 

(A.3) 

Here we show the energies measured from the one of the 
2a core calculated by superposing the same basis wave 
functions as the core in Eq. ([T]), that is, two as whose 
distances are 1 — 6 fm, and projected it to J T = + . In 
FIG. [HI it is shown that energies for most of states are al- 
most independent of the parameter S as shown by plateau 
lines in FIG. [6l These states may correspond to contin- 
uum states. Among those continuum states, there ex- 
ist resonance states whose energy decreases significantly 
as the parameter 5 increases. When 5 is large enough, 
the resonance states are decoupled from continuum states 
thoroughly to form bound states. In such a decoupling 



Appendix: Identification of the resonance states 



Here we show the additional way to identify the res- 
onance states in continuum spectra which we assign to 
the members of the K v = Q + bands labeled as 0^ 4 and 
so on in this paper. 



Since we superpose a number of wave functions to 
describe the quasi-bound cluster states, there appear 
a number of continuum states in addition to the res- 
onance states. As we describe in Sec. IIII1 we identify 
the 0j 4 states by means of the significant amplitude of 
a+a+dincutron and 6 He+a in a finite region respec- 
tively. However, these criteria may seem to be ambigu- 
ous. To make sure the identification further, we per- 
form additional analysis following the pseudo potential 



21|. 



method done in Ref. 

In this method, we introduce a pseudo potential, V, 




FIG. 6: The transition of the energies of the + states as 
the pseudo potential is increased. The horizontal axis is the 
parameter of the strength of the pseudo potential, 8. The 
points at 8 — correspond to the + states plotted in FIG.Q] 
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region, the 0^ and 0^ states are composed of the promi- 
nent components of ct+a+dineutron and 6 Hc+a, respec- 
tively. The corresponding states at 6 = (the points at 
S = in FIG. [5]) are consistent with the assignment of 
the 0^ and 0|" states discussed in Sec. Mil 

In the present calculation, main continuum states are 
2a core plus free two neutrons due to the present bound 
state approximation where 6 He+a continuum states are 
not considered sufficiently. 6 He+o; continuum states can 
be included by superposing 6 He+a cluster wave functions 
with larger d c . When ones up to d c = 15 fm are super- 
posed, some 6 He+a continuum states appear below the 
O4 state. We have applied the pseudo potential method 
with such bases and confirmed that the resonance states 
are decoupled from both continuum states of 2a+n+n 
and 6 He+a in the same manner in the above analysis. 

To assign the other members of their rotational bands, 
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